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Abstract 

We show that various identities from [1] and [3] involving Gould-Hopper polynomials 
can be deduced from the real but also complex orthogonal invariance of multivari- 
ate Gaussian distributions. We also deduce from this principle a useful stochastic 
representation for the inner product of two non-centered Gaussian vectors and two 
non-centered Gaussian matrices. 

Resume 

Formules de sommation pour les polynomes de Gould et Hopper Nous 
montrons que des identites pour les polynomes de Gould et Hopper issues de [1] et [3] 
peuvent etre deduites de l'invariance orthogonale reelle mais aussi complexe des lois 
Gaussiennes multivariees. Nous deduisons aussi de cette propriete une representation 
stochastique du produit scalaire de deux vecteurs gaussiens non-centres, ainsi que 
de deux matrices gaussiennes non-centrees. 
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1 An identity by Graczyck and Nowak 

We adopt the multi-index notation: with n > 1, a multi-index is denoted as 
m = (mi, • • • j m n) £ N n and its length |m| = mi + . . . + m n . In particular, the 
multi-index factorial is m! = mi! . . .m n \ and, with x e W 1 , the multi-index 
power x— = x™ 1 . . . x™ n . The Gould- Hopper polynomials [2j are 

Lm/2j , 

9m (x,p) = £ - ■ ,y x m ~ 2k ; x,peR. (1) 

Their multi-index version reads gm (x,£>) = IliLi fa • 
We use the following notations: X ~ A^(m y , Rx ) expresses the fact that the 
vector X is Gaussian with mean m x and covariance matrix Rx ', the expecta- 
tion operator is denoted as E, the Euclidean norm of vector z as \z\ and (n)- is 
the Pochhammer symbol T ^(n) > • Finally, underlined variables denote vectors, 
random variables are capitalized, x* denotes the transpose vector of x and ~ 
means equality in distribution. 

Our first result is a stochastic representation of the inner product of two non 
centered Gaussian vectors. 

Theorem 1.1 If x, G W 1 , pet and X, Y are independent random vectors 
~N(Q,L) then 

(x + y/pNf ( u + VpM) ~ (x + y/pNi) (y + JpMi) + pZ n _iN (2) 

where N, Ni and M 1 are independent standard Gaussian random variables, 

_ \x + y\ + \x~m\ _ \x + m\~ \x~u\ /on 
x- 2 ' y ~ 2 [) 

and Z n _i is chi- distributed with n — 1 degrees of freedom and independent of 
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N, Ni and M x . 

Proof 1.2 The proof is based on the polarization identity, that allows to ex- 
press this inner product as 



' \x + Vp~N+U+ VpM\ 2 - \x+ y/pN-y- y/pMi' 



\x + u+ \fzpK\ 2 -\x-u + \f^vM\^ 



where vectors M — (N+ M) and jV = ^ (N— M) are again Gaussian and 
independent as a consequence of the orthogonal invariance. 

The next step is again a consequence of this invariance: if £? ~ N{0, I n ) and 
z G W 1 , then the norm of the random vector z+ G depends on \z\ only; more 
precisely 



\z+G\ ~ | \4& + G\ = V(U| + Gi) 2 + G\ + . . . + G\ 

where e x is the first (or any) column vector of the identity matrix L^. We 
deduce that 

(x+ yfp0 (u+VpM) ~ t [lk + iz|ei + y/p~K+ y/pM\ 2 ~ \\x-uki + y/pR- y/p0 

and by the polarization identity, with X and Y as in Theorem this ex- 
pression simplifies to 

n 

(xe x + ^p-Nf (ye, + ^M) = (x + y/pNj (y + y/pM^ +pJ2N i M i . 

8=2 

The proof follows from the identity in distribution, with Z n _i ~ Xn-i indepen- 
dent of N ~ Af (0,1): 

n / n \ 2 

Y^NiM^ y£ M i) N = Z n . l N 

i=2 \i=2 ) 
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We deduce an elementary proof of the following [3] Thm 3], originally derived 
using generating functions. 



Theorem 1.3 For all M e N, p G R and x, y_ e IR n , we have, with x and y 
as in (Ej), 

1 LM/2J (2p) 2j (n-l\ 

E —,9m (x, p) g m (u, P) = E TTJ _ a\\ ~^T~ 9m-2 3 (x, p) g M -2j (y, p) (4) 
\m\=M — j=o J' \ Z JJ- v z 'j 

Proof 1.4 VFe use the moment representation of Gould-Hopper polynomials 
g m {x,p)=^[x + ^2pNy (5) 

where N ~ Af (0, 1): it is a consequence of definition (QJ) and of the fact that 
EN 2k = §§f and EiV 2fc+1 = 0. We t/ien recognize the left-hand side of (fjj as 
i/ie multinomial expansion 



|TO|=M — i=l v 7 v 

= ^ E (E + v^*) (w + V^)) 



where M and N ~ JV f 0, J n J fl^e independent. The sum is identified as the 
inner product 

Using the multinomial expansion of the M — th power of the stochastic equiv- 
alent of this expression and taking expectations with ~EZ 2 _i = 2 J (j^^j > 
we obtain the desired result. 
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2 Two identities by Daboul and Mizrahi 

The group of complex rotations Oq (n) is the set of n x n complex matrices Q_ 
such that £2£2' = Qt (2 = [n- The real inner product x* y_ = J27=i x iViy ^ X 
C n , is preserved under the action of this group. In p], sum rules for Gould- 
Hopper polynomials are proved using the covariant transformation property 
under Oc (n) of the raising operators associated with these polynomials. We 
show here that these sum rules can be equivalently deduced from the following 
complex rotational invariance principle. 

Theorem 2.1 For any Q E Oc(n) and N ~ A/^O,/^), the Gould-Hopper 
polynomial of degree m reads 



9m(x,p)=E\x + J2p[QN) ) , Vje[l,n 



j , 



Proof 2.2 We only need to check that y(2NJ ~ A/"(0, 1); but, using the fact 
that the identity Eexp (2ttNj) = exp (— t 2 ) holds for any complex value oft, 
the characteristic function of this random variable reads 

Eexp [it (QNj J = exp (^-t 2 £ 0| fc J = exp (-t 2 ) 
by the complex orthogonality condition. 

Using this probabilistic invariance, we derive a simple proof of the summation 
theorem [1, Prop. 1] 

Theorem 2.3 For any Q E Oc (n) , 

g m (((2x),p)= Yl ( m J fi^dmj (xj,p) 

Proof 2.4 The proof is based on the complex rotational invariance (Theorem 
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\2.1\) as follows 



9m ((Qx) ,p)=E ((ax) + ^2p (qn) y ~ E [j^Oij (xj + ^TpN j 



\j=l 

The desired result is obtained by expanding this power using the multinomial 
formula as 

e e L m m )fi°7U+M^r = e L m Jn°?^M 



|m|=m \ mi ' • • • ' m? v 3=1 v y \m\=m V m i' • • • ' m ™/ j=l 



We now prove equally easily the following generalized factorization sum rule 
PQ Prop.3] 

Theorem 2.5 With p e E, c, s, x, y G C and c 2 + s 2 = 1 ; 

mi+m2 

g mi (cx- sy,p)g mi (sx + cy,p) = ^ Cr^^.r (c, s) # r (x,p) 9m,i+m,2—r (y,p) 

r=0 



where coefficients C mitTn2tr (c, s) are given by 

C„,r (C, 8) = h) (~ir- r+l <?»+r-" 8 ™-++* ( 6 ) 



=0 



Proof 2.6 By the complex orthogonal invariance, [JVi, iVa] ~ A/"(0, ^2) 
[ciVi - siV 2 , siVi + ciV 2 ] so that 



g mi (cx - sy,p) g m2 (sx + cy,p) 

= E I cx - sy + sj2p (ciVi - sN 2 ) J E (sx + cy + y2p (sNi + cN 2 ) J 

= E E E (7) (?) (- 1 ) mi " fc c m2+fc -^ mi+/ ~ fc (x + v^iVi)^ (y + vW m "' +n_fe 

mi+mi 

= C miim2 ^ r (c, s) g r (x,p) g mi + m2 - r (y,p) 

with coefficients C miim2jr (c, s) as in (0|) ; which coincides with £ZJ eg. 
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We note that all former sum rules apply to the Hermite polynomials since 
H n (x) = g n \ 2x, a/— 2^ . Moreover, theorem 11.11 extends to the matrix case as 
follows, with ||20 f denoting the Frobenius norm: 

Theorem 2.7 If z,y_ G R mxn and M, N e K mxn are independent matrices 
with i.i.d. standard complex Gaussian entries then, with Mi, Ni, N ~ J\f (0, 1) 
and Z mn _i ~ Xmn-i independent, 

tr(g + Kf{g + M) ~ (x + iVi) (y + Mi) + Z ron _iiV 

wt/i x = l/2(||x + ^|| F + ||s -g\\p),y= l/2(||g + y||| - lis - ^|||) ■ 

Proof 2.8 Denoting the i—th column vector of z and u,v the vectors with 
components Ui, Uj = l/2(||xj+^j|| 2 d=||3: i — ^||), we deduce that 1/2 (||« + j;|| 2 ± \\u — v\\) = 
x,y and by theorem LPl 

m m 

tr(z + Kfiu + M) = Efe + ^)'G& + ~ £ + ^) + M *) + xS-i^. 

Using again theorem the first sum is distributed as (x + Ni) (y + Mi) + 
Xm-iG and since + Y%Li Xn-iQ ~ Xnm-iN, the result follows. 
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